Competitively Coupled Maps and Spatial Pattern Formation 
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Spatial pattern formation is a key feature of many natural systems in physics, chemistry and 
biology. The essential theoretical issue in understanding pattern formation is to explain how a 
spatially homogeneous initial state can undergo spontaneous symmetry breaking leading to a stable 
spatial pattern. This problem is most commonly studied using partial differential equations to model 
a reaction-diffusion system of the type introduced by Turing. We report here on a much simpler 
and more robust model of spatial pattern formation, which is formulated as a novel type of coupled 
map lattice. In our model, the local site dynamics are coupled through a competitive, rather than 
diffusive, interaction. Depending only on the strength of the interaction, this competitive coupling 
results in spontaneous symmetry breaking of a homogeneous initial configuration and the formation 
of stable spatial patterns. This mechanism is very robust and produces stable pattern formation for 
a wide variety of spatial geometries, even when the local site dynamics is trivial. 

PACS numbers: 05.45.-a, 89.75.Kd, 89.75.Fb 



I. INTRODUCTION 

Achieving a satisfactory understanding of spatial pat- 
tern formation is an enduring problem in the physical and 
biological sciences [Ij]. At the fundamental level the es- 
sential issue is to explain how a spatially inhomogeneous 
stationary pattern can emerge dynamically from a spa- 
tially homogeneous initial state. The most widely stud- 
ied mechanism for producing such spatial patterns is that 
proposed by Turing Turing's model is formulated as 
a system of coupled reaction-diffusion equations that de- 
scribe the dynamical interaction of two chemical species 
or morphogens. Pattern formation in Turing's system 
results from the counterintuitive fact that, under suit- 
able conditions, the diffusion of the morphogens can drive 
symmetry breaking in the initial homogeneous configura- 
tion. This is possible in Turing's system because diffusion 
is acting in concert with other processes: namely, one of 
the morphogens (which undergoes an autocatalytic re- 
action) activates the other morphogen, while this mor- 
phogen in turn inhibits the former. Whether such a 
system can produce stationary spatial patterns (as op- 
posed to, for example, spiral waves as in the Belosov- 
Zhabotinsky reaction [3]) depends delicately on the dif- 
fusion rates of the morphogens, and on the rate at which 
the inhibitor responds to changes in the activator con- 
centration y . This subtle dependence of the mechanism 
on parameter values makes Turing structures challeng- 
ing to produce experimentally [5], which suggests that 
the mechanism may have limited applicability to pattern 
formation in many real- world situations. 

A variety of spatio-temporal patterns have also been 
shown to arise in the context of coupled-map lattices. 
Originally introduced by Kaneko, these well-studied sys- 
tems [ol consider a lattice geometry where the temporal 
site dynamics of a single variable is described by a dis- 
crete time mapping, while a functional dependence on 
the neighboring values of the site variable determines the 
spatial dynamics. This coupling between sites is often 



diffusive in character [6|, though a number of variants 
have been studied (see for example Ref. [7]). A com- 
mon feature of these models is that the uncoupled site 
dynamics needs to be sufficiently complex (often chaotic 
in terms of the local map dynamics) in order to generate 
non-trivial spatial patterns on introducing spatial inter- 
actions. As is clear from the case of diffusive coupling, 
the spatial coupling serves to coarse-grain the large varia- 
tions at the sites leading to patterns. A natural question 
that proceeds from this is whether or not there are mech- 
anisms to generate non-trivial spatial patterns when the 
on-site dynamics is fully trivial (i.e. the local dynamics 
has a stable equilibrium point)? 

The purpose of this paper is to study a model of spatial 
pattern formation which is both simpler and more robust 
than Turing-type models and readily produces station- 
ary patterns from a homogeneous initial configuration. 
Our model is formulated as a variant of a coupled map 
lattice [6], in which there is a single dynamical variable 
(i.e. a single species or morphogen level) at each spa- 
tial site and the site dynamics are coupled through a 
competitive (that is, inhibitory) interaction rather than 
through a conventional diffusive coupling. The competi- 
tive coupling in our model allows the production of com- 
plex spatial patterns even when the dynamics at each 
site is trivial, in the sense that the local dynamics ex- 
hibits a stable fixed point. This is in stark contrast to 
the dynamical behavior in conventional diffusively cou- 
pled map lattices where trivial site dynamics can only 
result in a spatially homogeneous state [8-11]. This lat- 
ter result reflects the fact that for a single species model 
diffusion can never have a destabilizing effect. In diffu- 
sively coupled map lattices, pattern formation can occur 
but this is only possible when the local site dynamics is 
complex. Indeed, non-trivial pattern formation typically 
requires that the local dynamics is chaotic, [g] 

In order to appreciate the importance of competitive 
coupling, as opposed to diffusive coupling, for spatial pat- 
tern formation in map lattices, it is necessary to briefly 
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review the stability properties of diffusively coupled map 
lattices. An important issue in formulating a diffusively 
coupled map lattice is the form taken for the diffusive 
coupling term. The initial work on coupled map lat- 
tices [g^, 8] considered maps of the form 

Xi{t^l) = f{Xi{t))^d[x^^l{t)^Xi-^{t)-2Xi{t)], (1) 

in one-dimension with nearest neighbor interactions, and 
the obvious generalizations to other spatial dimensions 
and interaction neighborhoods. The map f{x) defines 
the local dynamics, and d represents a "diffusion" con- 
stant. The case in which f{x) is the logistic map 
f{x) = rx(l — x) has often been considered in the litera- 
ture [H, m, [ill . The question of whether non-trivial spatial 
patterns can arise when the local dynamics is trivial (i.e., 
when the local dynamics f{x) has a stable fixed point) 
depends on the stability of the coupled map ([T]). For 
f{x) of the form of the logistic map the stability theory 
of ^ has been derived in Ref. [8], and also in Ref.[10]. 
It has been shown that ^ always has a stable homoge- 
neous fixed point if the map f{x) has a stable fixed point. 
Thus, if the local dynamics is trivial then it is impossi- 
ble to obtain non-trivial spatial pattern formation in the 
coupled map lattice. For this reason attempts to find 
non-trivial patterns in coupled maps of this form have 
used local maps with non-trivial dynamics. For example, 
in Ref. [l2| spatial pattern formation was found in the 
two-dimensional version of ([T]) with f{x) = rx(l — x), 
and r = 3.25 and r = 3.30, which correspond in both 
cases to the local dynamics being a stable 2-cycle. 

The coupled map lattice ^ has, however, a funda- 
mental problem, which was first identified in [Tsl (see 
also [1Q|)' the coupled map ([1]) is mathematically not 
well-defined — namely, the iterates of ([1]) can become 
unbounded for perfectly reasonable values of d and ini- 
tial values Xi{0). This lack of boundedness is an indi- 
cation that a coupled map lattice of the form ([1]) does 
not provide an appropriate model of a discrete reaction- 
diffusion system. In fact, a careful analysis shows that 
a coupled map of the form ([1]), for any local dynamics, 
spatial dimension, and interaction neighborhood, does 
not have a coherent physical interpretation in terms of 
a discrete analogue of a reaction-diffusion process prj . 
Systems governed by coupled maps of the form ([T]) have 
the unphysical property that individuals (e.g., molecules 
or organisms) may disintegrate (or die) and still diffuse 
(or disperse) (see for a detailed discussion). This 

lack of mathematical well-definedness together with the 
lack of any coherent physical interpretation has resulted 
in coupled map lattices with diffusive interactions of the 
form taken in ([!]) falling out of use. 

The solution to the problems alluded to above is simple 
and elegant — to separate the reaction (or reproduction) 
process represented by the local map f{x) from the diffu- 
sive process [10, 11]. Each iterate of the coupled map is 
now viewed as having two stages. In stage 1 (the reaction 
stage) we have 

x',{t) = f{xi{t)). (2) 



This is followed by stage 2 (the diffusion stage) to give 
the final value of the next iterate 

x,(t + 1) = x[{t)^d[x[^^{t)^x[_^{t) - 2x[{t)]. (3) 

Combining these two stages we obtain the complete cou- 
pled map lattice 

Xi{t+1) = f{xiit))+d[f{xi+i{t))+fixi.i{t))-2f{xim- 

(4) 

This form of the coupled map lattice generalizes directly 
to arbitrary spatial dimensions and interaction neighbor- 
hoods. This form (j4]) of the coupled map lattice com- 
pletely resolves the two problems associated with the ear- 
lier form of the coupled map lattice ([1]) : it is easy to show 
that iterates of dH are always bounded for reasonable 
values of d and ^^(0) [10], and (j4]) has a clear physical in- 
terpretation as a discrete reaction-diffusion system with 
a reaction stage followed by a diffusion stage [10]. 

Due to these advantages all recent work on diffusively 
coupled map lattices have used the form dU , and its gen- 
eralizations. The stability theory of coupled map lattices 
of the form @ has been studied in great generality. It 
has been shown that for such a map lattice with arbi- 
trary local map /(x), spatial dimension, and interaction 
neighborhood, if f{x) has a stable fixed point then the 
corresponding coupled map has a stable homogeneous 
fixed point [sllili- Thus, if the local dynamics is trivial 
it is impossible for any coupled map lattice of the general 
form of (H)) to produce non-trivial spatial patterns. We 
note that this general result is in accord with physical 
intuition: if the local map approaches the same equilib- 
rium value at every site, then it would be most surprising 
if the introduction of a diffusive spatial coupling (which 
will act to reduce differences between site variables at 
different positions) could result in a spatially inhomoge- 
neous equilibrium. It follows from these stability results 
that non-trivial global dynamics is only possible in any 
coupled lattice map of the form of dH if the local dynam- 
ics defined by f{x) is non-trivial, in the sense that it does 
not have a stable fixed point. It is for this reason that 
any coupled map of the form @ that displays non-trivial 
global dynamics is necessarily based on a local map that 
has non-trivial dynamics. Examples of such situations 
include: Ref. [l^, in which the local dynamics is defined 
by a chaotic tent map; Ref. [15], where the local dy- 
namics is given by a chaotic logistic map; and Ref. [i6| 
in which the local map has a stable 3-cycle. 

The notion of coupling local site dynamics via a com- 
petitive interaction, as in our model, has naturally arisen 
in both theoretical ecology and developmental biology. In 
the ecological context, the dynamical variable Xi{t) can 
be interpreted as the size of the population at site i and 
at time t. Since the growth in the population at a given 
site will be limited both by competition with individuals 
at that site and also with individuals from neighboring 
populations, it is natural to consider competitive cou- 
pling in the population dynamics of such systems. For 
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example, Hara and co-workers [17] addressed the effects 
of local versus global competition on size distributions in 
tree populations. The site dynamics there was described 
by a single logistic differential equation, which exhibits 
only a stable equilibrium point, which corresponds to a 
uniform tree size. However, on including either global 
or local spatial competition, there emerges a parameter 
threshold beyond which variation in size distribution is 
seen. A discrete-time model of spatial population dy- 
namics with local competition has also been addressed in 
Ref. [18]. The notion of competitive coupling as a basic 
mechanism for pattern formation is also closely related 
to the idea of lateral inhibition, in which the growth of a 
structure at a given location inhibits the growth of sim- 
ilar structures in a region surrounding the focal object, 
which has been studied in connection to pattern forma- 
tion in developmental biology [19]. 

It is interesting to note that a similar model to the 
one we study here has been investigated in Refs. (2Q| - 
[22[ in the context of understanding income distributions. 
Although the model studied in these papers is formally 
similar to our model, the aims of our work and that of 
Refs. [2Q|-[22| are very different. Whereas our focus is 
on how a homogeneous initial state can undergo spon- 
taneous symmetry breaking and result in stable pattern 
formation, Refs. |2 01-12 21] do not contain any discussion 
of spatial pattern formation and are, instead, concerned 
with modeling wealth or income distributions in society. 

It is also worth mentioning that the form of our model 
is typical of that of various discrete competition mod- 
els jlsl, [23 . The novel aspect of our work is that we show 
that a simple, coupled map model incorporating spatial 
competition can produce stable spatial pattern forma- 
tion, even when the local dynamics is completely trivial. 
In contrast, the models studied in [23|,[23 are not spatial 
models and are unconcerned with any aspects of spatial 
pattern formation. 

Finally, in contrast to the specific (and in some cases 
rather complicated) models of spatio-temporal dynamics 
with local competition that have been previously investi- 
gated, our aim in this paper is to explore in some detail 
a very simple class of discrete-time spatio-temporal dy- 
namical models in which competitive interactions lead 
to spatial pattern formation. We show that in this gen- 
eral class of models, pattern formation is a robust conse- 
quence of local spatial competition, and does not depend 
on the specifics of either the local dynamics or the spatial 
geometry. 



II. COMPETITIVELY COUPLED MAPS 

We begin, for the sake of generality, by defining a com- 
petitively coupled system of maps for an arbitrary net- 
work, which we shall refer to as a competitively coupled 
map network (CCMN). However, most of this paper will 
focus on the case of relevance for pattern formation in 
which the network is a lattice, corresponding to compet- 



itively coupled map lattices (CCML). 

We define a CCMN as follows. First, let / : R ^ R be a 
smooth, one-dimensional, unimodal map. Further, let us 
assume that / can be written in the form f{x) = xF{x) 
where F : R ^ R is also a smooth map that is monotoni- 
cally decreasing on the domain of interest. We note that 
many celebrated maps / can be expressed in this form. 
Well-known examples [25] include the logistic map where 
F{x) = r(l— x), the Ricker map where F{x) = Xe~^^ ^ the 
Hassell map where F{x) = A/(l + ax)^, and the Maynard 
Smith map for which = A/(l + ax^). The parame- 
ters r. A, a and b are positive constants. The parameter 
a can be scaled away so we can, without loss of general- 
ity, set a = 1. It is well-known that all these maps have 
qualitatively identical dynamical features, including sta- 
ble fixed points, limit cycles and a period-doubling route 
to chaos [25| . Here we focus on pattern formation ariseing 
from the Ricker, Hassell and Maynard Smith maps, which 
are a representative sample of one-dimensional, unimodal 
maps for which the iterates are always non- negative. We 
do not consider the logistic map as the possibility of the 
iterates becoming negative results in technical complica- 
tions which have no bearing on the pattern formation 
mechanism we are investigating. Fut her more, as we wish 
to emphasize the role competitive interactions play in 
pattern formation, we chose to employ maps with the 
simplest possible dynamical behavior. As such, we al- 
ways consider here only parameter values for which the 
maps have a stable fixed point. We reiterate that for such 
parameter values, diffusively coupled map lattices always 
evolve to a spatially homogeneous state and pattern for- 
mation is impossible. Thus, any pattern formation that 
results in our model is a consequence of local competition 
and not the complexity of the local dynamics. 

First we consider in some detail the case in which the 
local dynamics is given by the Ricker map and then we 
will expand the discussion to include the Hassell and 
Maynard Smith maps. For the Ricker map, f{x) = 
Axe~^, the fixed points are x = and = InA given 
by the condition = 1. The corresponding stability 

is determined by the function (j){x) = F{x) + xF'{x). A 
fixed point x* is linearly stable if |0(x*)| < 1 and unsta- 
ble if |0(x*)| > 1. For the Ricker map, x = is stable 
for < A < 1 while is stable for |1 — In A| < 1, that is 
for 1 < A < e^. When A exceeds 1, the stable fixed point 
X — loses stability through a transcritical bifurcation 
and gives birth to the stable fixed point = In A. Here 
we always restrict attention to A < corresponding to 
stable fixed point dynamics. 

To form a CCMN based on a map f{x) = xF{x)^ we 
first introduce a network F, i.e. F is a simple, undirected 
graph with m vertices (or nodes) which are labeled by 
i = 1, • • • , m. The topology of F is determined by the ad- 
jacency matrix A = (aij)^ i, j = 1, • • • , m, where aij = 1 
if nodes i and j are connected by an edge (or hnk) and 
aij = otherwise. Given a node z G F we define the set 
of neighbors of i to he N{i) = {j G F : a^j = 1}. The 
CCMN associated with / and F is a discrete-time dynam- 
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ical system ^ on W^. The state of this system at time t is 
determined by a vector x(t) = • • • G M"^, 

where {t) represents the value of the state variable asso- 
ciated to node i of the network at time t. The dynamical 
system ^ gives the time evolution of the state x(t), i.e. 
^(x(t)) = x(t + 1), where <l> is defined by: 



j e N{i) 



(5) 



for all z G r. In this expression a is a non- negative 
parameter that represents the strength of the competitive 
interaction between site i and a neighboring site j. 

The first step in understanding the dynamics of (jSj) is 
to study its fixed points and their stability. In discussing 
the fixed points of $ it is convenient to assume that the 
network F is a regular graph of degree /c, which means 
that every node has exactly k edges incident upon it. In 
this case, there are two homogeneous fixed points of ([5]), 
where this term refers to every node z G F having the 
same value for the state variable. The first of these is 
the trivial fixed point x = = (0, 0, • • • ,0) while the 
non-trivial fixed point is x* = • • • where 

satisfies F[{1 + ka)x*] = 1. For the case of the Ricker 
map, this latter fixed point is given by x* 



In A 
l-\-ka ' 



The stability of the fixed points is determined by the 
eigenvalues of the Jacobian matrix J(x) which has ele- 
ments Jij (x) = ^^^^^^ . Direct calculation from (|5|) shows 
that the Jacobian evaluated at the trivial fixed point 
X = is J(0) = F{0) ■ /, where / is the m x m iden- 
tity matrix. Thus, every eigenvalue is equal to F(0), 
which means simply that the eigenvalue of the original 
uncoupled map determines the stability at the fixed point 
X = 0. Therefore, the CCMN dynamics is trivial if and 
only if the dynamics of the uncoupled map / is stable for 
X = 0, which for the Ricker map means A < 1. 

Considering now the non-trivial, homogeneous fixed 
point X*, we find that the Jacobian evaluated at x* is 



(6) 

where A is the (m x m) adjacency matrix of F. The 
stability of x* is determined by the m eigenvalues ipj of 
J. 

It follows immediately from eq.® that J(x*) is a circu- 
lant matrix whenever A is a circulant. Note that circulant 
matrices are fully specified by the first row vector as all 
other rows are cyclic permutations of the first row where 
the offset corresponds to the row index [26]. Those net- 
works for which the adjacency matrix is a circulant form 
an interesting class, known as circulant graphs. The net- 
works we study here are all either circulant graphs or are 
networks that can be directly related to circulant graphs. 
For such networks, the eigenvalues of J(x*) can be di- 
rectly computed using (|6|). 



For the Ricker map, the non-trivial fixed point is x* = 
^ and the network Jacobian (j6|) evaluated at this fixed 



point is 



J(x* 



In A 



■I 



a In A 

ka-\-l 



■A. 



(7) 



Given any circulant graph F it is possible to explicitly 
compute the eigenvalues V^j, j = 0, ••• ,m — 1, of this 
Jacobian from which the stability of the non-trivial fixed 
point in the CCMN can be established. 



III. SPATIAL PATTERN FORMATION 

Within the general framework of a CCMN, we will now 
focus on such systems in the cases of most interest for 
pattern formation, when the network F is a spatial lat- 
tice in one or two dimensions. We begin with the simpler 
case where F is a one-dimensional lattice with m nodes 
subject to periodic boundary conditions. Let us also re- 
strict the analysis to nearest neighbor interactions where 
site i interacts with sites — 1) and (i + 1). The gener- 
alization to other interactions is straightforward. Thus, 
in this case, k = 2 and the homogeneous fixed point 
of the competitively coupled map lattice is 



In A 



2a+l • 

Given the periodic boundary conditions, the Jacobian (|6]) 
is a symmetric circulant matrix and the associated real 
eigenvalues can be explicitly computed from the prop- 
erties of circulant matrices [26]. This calculation gives 
ipo = \nX = (/>(a:*), which is the eigenvalue of the uncou- 
pled Ricker map evaluated at the non-trivial fixed point 
X* = InA of the map. Thus, |V^o| < 1 whenever the fixed 
point X* = In A of the Ricker map is stable. However, if 
X* = In A is unstable for the uncoupled Ricker map, then 
l^/^ol > 1 and the homogeneous fixed point x* = ^^^^ of 
the CCML is also unstable. 

Consequently, stability of the fixed point = In A is a 
necessary condition for the stability of the homogeneous 
fixed point x* of the CCML, but it is not sufficient. This 
may be seen by computing il^rn/2i foi" ^ even, or '0(m-i)/2 7 



for odd m. For even m we obtain '0^/2 = 1 



2a-l 



InA, 



2cx+l 

which is the same result as that for the odd m case, 
'0(m-i)/2 7 for m >> 1. It is straightforward to show 
from the eigenvalues of J(x*) that if |V^o| < 1, (i-e. for 
1 < A < e^), then the stability of the homogeneous fixed 
point X* is determined by = 1 + In A, for m >> 1. 

Thus, X* is stable if < 1, i.e. for a < 1/2 and un- 
stable for a > 1/2. This clearly demonstrates that the 
competitive coupling results in the stability of the non- 
trivial homogeneous fixed point being distinct from the 
stability of the nodal Ricker map dynamics. Therefore, for 
1 < A < (i.e. for trivial site dynamics), a > 1/2 re- 
sults in spontaneous symmmetry breaking of the spatially 
homogeneous initial state. As we shall show, this results 
in the emergence of complicated stable spatial patterns. 
In this one dimensional situation, the dynamics at a 
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site i given by (j5j) takes the form 

Xi{t + 1) = Xx.(^t)e-^Mt)+c.i^^-lit)+x.^ 



lit))] 



(8) 



It is apparent from (|8]) that our CCML naturally in- 
corporates the principle of local activation and longer- 
ranged inhibition that is the essential, but indirectly re- 
alized, mechanism of pattern formation in Turing-type 
systems [4]. It is also clear from the negative exponen- 
tial in dHJ that a site surrounded by neighbors with large 
values of the site variable x will feel strongly the effects 
of the coupling, and will be suppressed. However, while 
a strong coupling has an inhibitory effect on a site be- 
tween two neighbors with large x values, the same cou- 
pling can leave the larger neighbors relatively unaffected, 
since low surrounding sites have little inhibitory effect 
on their neighbors. This is in strong contrast to diffusive 
coupling [6], which for a single species can only act to 
uniformize the levels of the site variable [§1,0. It is this 
difference which is responsible for the pattern formation 
seen in our model. 

Figure [1] illustrates the change in stability in our model 
at a = 1/2 for a one-dimensional chain with periodic 
boundary conditions. In Fig. l(a)[ the differences in 
initial populations at the sites are quickly erased and 
the stable, non-trivial, spatially homogeneous state is 
reached. By contrast, on crossing the threshold a = 1/2, 
a ve ry in teresting spatial pattern emerges as seen in 
Fig. l(b)[ We note that at early times, the alternat- 
ing pattern discussed earlier is clearly visible which then 
evolves to a more complex pattern in time. This tran- 
sition happens progressively later in time as one ap- 
proaches the a = 1/2 threshold. It should be emphasized 
that for this value of the map parameter A = 6 < e^, the 
uncoupled site dynamics would converge to a fixed point 
which clearly demonstrates that the complexity seen in 
the spatial pattern is an effect of symmetry breaking. 

Let us now consider F to be a two-dimensional square 
lattice with nearest neighbor interactions, where each site 
(i, j) is coupled with sites (z — 1, j), (z, j — 1), (i, j + 1) and 
(i + 1, j). The map dynamics is given by the appropriate 
specialization of ([5]). Our earlier analysis can be readily 
extended to this situation to show that a = 1/4 is the 
critical value of the inter-site coupling. 

For A < e^, the Ricker map evolves to a stable fixed 
point, X* = In A. For these A values and with coupling 
a < 1/4, the homogeneous fixed point x* of the CCML 
is stable and the value at each site converges to In A/(l + 
4a). When the coupling reaches a = 1/4 we see the 
first qualitative change in the dynamics of the CCML, 
a "checkerboard" pattern consisting of alternating high 
and low values. The checkerboard pattern arises as a 
consequence of the short range activation and long-range 
inhibition inherent in competitive coupling. 

At a = 1 /4 the differences in neighboring sites evolves 
slowly, and the checkerboard emerges after several tens 
of thousands of iterations. However, with just slightly 
larger coupling {a = 0.26), the formation of the checker- 
board becomes much more rapid. In particular, the 
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FIG. 1. Evolution of the population with time t on a one- 
dimensional chain with 200 nodes. In (a), the coupling 
a = 0.49 which is smaller than the threshold value. Panel 
(b) corresponds to a = 0.53 and results in a non-trivial pat- 
tern. In both cases, the parameter A = 6 for which the site 
dynamics has a stable fixed point. 



checkerboard pattern forms at several different regions 
of the lattice at the same time. This leads to the forma- 
tion of "domain walls" between neighboring checkerboard 
patterns when the high and low sites are out of phase. As 
seen in Figs. |2(a)[ the boundaries of these domains either 
wrap around the toroidal geometry or form simple closed 
loops on the lattice. At these values of the coupling, 
simple loops on the lattice are transient patterns which 
ultimately shrink to a point, over times which are larger 
than 10^ iterations. By contrast, the toroidal patterns 
are stable and may be viewed as topological defects as- 
sociated with the non-trivial first homotopy group of the 
torus. As a increases, the minimum and maximum val- 
ues of the underlying checkerboard pattern do not vary 
significantly. However, the increased coupling strength 
results in the checkerboard emerging more rapidly and 
at more points on the lattice, resulting in complicated 
domain wall structures. Additionally, with the stronger 
coupHng, isolated closed loops no longer decay but b e- 
come as persistent as the toroidal patterns. (Fig. 2(b)). 



Having illustrated the ability to generate stable spa- 
tial patterns, we now include the effects of next-nearest 
neighbor interactions, where site (i^j) interacts with all 
eight sites surrounding it. The threshold value of a can 
be verified to remain at 1/4. We observe that in the next- 
nearest neighbor CCML the checkerboard pattern is no 
longer a feature of the model. While the local activation 



(a) (b) 
FIG. 2. Typical patterns seen in two-dimensional, nearest 
neighbor CCMLs, illustrated for the Ricker map. In both 
cases, A = 5, = 100^, (a) a — 0.26 and shows both toroidal 
patterns and a long-term transient loop, (b) a — 0.35 
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FIG. 4. Typical patterns in next-nearest-neighbor CCMLs. 
A = 5 iV = 100^ and (a) a = 0.26 and (b) a - 0.35. 



and long range inhibition still persist, next nearest neigh- 
bor interactions allow or a new possibility, namely, the 
occurence of geometrical frustration [27]. The resulting 
structure closely resembles glassy spin systems [28^ . 
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FIG. 3. Schematics illustrating (a) the basic checkerboard 
pattern arising with nearest-neighbor interactions and (b) the 
disappearance of the effect due to geometrical frustration once 
next-nearest neighbor interactions are included. The + and 
— symbols refer to high and low site variable values, respec- 
tively. The symbol ± in the central cell of (b) indicates that, 
with next-nearest neighbor interactions, it is not possible for 
all nine site variables in the fundamental 3x3 neighborhood 
to assume their dynamically preferred values. This obstruc- 
tion is the origin of the geometrical frustration that occurs in 
CCMLs with next-nearest neighbor interactions. 

Local competitive interactions tend to result in sites 
with high values of the dynamical variable having neigh- 
bours with low values, and vice versa. With nearest- 
neighbor interactions there is no geometric obstruction 
to this occur ing as shown in Fig. |3(a)[ However, next- 
nearest-neighbor interactions introduce a geometrical ob- 
struction to such a dis tribution of dynamical site values, 
as shown in Fig. |3(b)[ The occurrence of this geometri- 
cal frustration results in the appearance of domains with 
distinct orientation as shown in Fig |H As seen from the 
second of the two panels shown, the domains get smaller 
with increasing a. We note that in all the cases shown, 
qualitative aspects of the results are independent of the 
initial conditions. Specifically, even very small (one part 
in 10^) initial spatial variations evolve dynamically into 
the spatial patterns shown. 

The pattern formation mechanism based on competi- 
tively coupled maps is applicable to any spatial geometry. 



In Fig. Owe show typical spatial petterns that result on 
triangular and hexagonal two-dimensional lattices, with 
nearest neighbor interactions and local Ricker map dy- 
namics. The pattern formation that occurs in these lat- 
tices is quantitatively very similar to that which takes 
place on square lattices with nearest-neighbor interac- 
tions. This illustrates the important result that pattern 
formation in our system results from the competitive in- 
teraction and is not dependent on the details of the spa- 
tial geometry. 





(a) 



(b) 



FIG. 5. Typical patterns in nearest-neighbor CCMLs on (a) 
hexagonal and (b) triangular lattices. A = 5, = 100^ and 
a — 0.45 in both cases. In each of these situations, the critical 
value of a is 1/3. 

Most of our discussion of pattern formation thus far 
has focused on local dynamics defined by the Ricker 
map. However, as we now show, pattern formation oc- 
curs in CCML when the local dynamics is described 
by any unimodal map. We illustrate this by consider- 
ing local dynamics defined either by the Hassell map, 
f{x) = \x/{l + ax)\oT the Maynard Smith map, 
f{x) — \x/{l + ax^) [25]. Here again, we restrict our- 
selves to parameter values for which these maps have only 
stable fixed points so as to make clear that any pattern 
formation is an outcome of local, competitive, spatial in- 
teraction and not a result of any dynamic complexity of 
the local map. 

The stability analysis of the CCML with Hassell or 
Maynard Smith maps is fully analogous to the Ricker 
map case discussed. Thus, it may be shown that, for 
a one-dimensional lattice with nearest-neighbor interac- 
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(a) (b) 
FIG. 6. Evolution of the population with time t on a one-dimensional periodic chain with 200 nodes. In (a), the Hassell map 
governs the local dynamics with b — 2 and A = 9 while panel (b) uses the Maynard Smith map for parameter values of 6 = 2 
and A = 4. The inter-site coupling a — 0.53 in both cases which is the same value used for the Ricker map case in Fig. [T]and 
it results in similar non- trivial patterns. 



tions, the critical value of a remains a = 1/2. Similarly, 
the critical value for both Hassell and Maynard Smith 
maps on a two-dimensional square lattice with nearest- 
neighbor interactions is a = 1/4. We show in Fig. [6] 
pattern formation for CCML with Hassell and Maynard 
Smith local dynamics in a one-dimensional latice geom- 
etry. The results are qualitatively similar to those seen 
earlier with Ricker map local dynamics. 




(a) (b) 

FIG. 7. Typical patterns in nearest-neighbor CCMLs with 
(a) Hassell and (b) Maynard-Smith local map dynamics. The 
relevant parameter values are (a) 6 = 2 and A = 9 and (b) 
5 = 2 and A = 4. A square lattice with N — 100^ sites and 
a = 0.26 was considered in both cases. 

In Fig. [71 pattern formation for the case of a square 
lattice with Hassell and Maynard Smith local dynamics 
are shown with nearest-neighbor interactions, which are 
again very similar to those seen earlier for the Ricker 



map dynamics. Adding next nearest neighbors results 
once again in the appearance of geometric frustration 
in the dynamics. These results make clear that pattern 
formation in CCML is a robust property which results 
from local competitive interactions and not the specifics 
of either the spatial geometry or the local dynamics. 
IV. CONCLUSIONS 

We have studied here a new model of spatial pat- 
tern formation based on a novel type of map lattice, in 
which the fundamental spatial interactions are competi- 
tive, rather than the diffusive ones that are typically con- 
sidered. This model directly incorporates the principle of 
local activation and long-range inhibition which is funda- 
mental to pattern formation in Turing systems, but which 
is only indirectly implemented via reaction-diffusion in- 
teractions in such systems. We have demonstrated that 
competitive coupling of spatially distinct maps results in 
complex stable patterns even when the site dynamics is 
trivial, in the sense that the local map has a stable fixed 
point. Further, we have shown that the patterns formed 
by this mechanism are qualitatively independent of the 
detailed form of the local site dynamics and of the spatial 
geometry. Thus, competitively coupled maps provide a 
simple and robust model of spatial pattern formation. 

B.S. is grateful to Rainer Scharf and Permanand Indie 
for helpful comments. 



[1] A. Gierer and H. Meinhardt, Kybernetik 12, 30 (1972); 
S. Douady and Y. Couder, Phys. Rev. Lett. 68, 2098 
(1992); A. Koch and H. Meinhardt, Rev. Mod. Phys. 66, 
1481 (1994); E. Boissonade, E. Dulos and P. De Kepper, 
in Chemical Waves and Patterns^ eds. R. Kapral and K. 
Showalter, 221 (Kluwer, Dordrecht, 1995). 

[2] A. Turing, Phil. Trans. R. Soc. B237, 37 (1952). 



[3] R. Kapral and K. Showalter (eds.), Chemical Waves 
and Patterns (Kluwer Academic Publishers, Dordrecht, 
1995); S. Jakubith, H. H. Rotermund, W. Engel, A. von 
Oertzen and G. Ertl, Phys. Rev. Lett. 65, 3013 (1990); G. 
Ertt, Science 254, 1750 (1991); J. Lechleiter, S. Girard, 
E. Peralta and D. Chapman, Science 252, 123 (1991); K. 
Agladze, J. P. Keener, S. C. Muller and A. Panfilov, Sci- 



8 



ence 264, 1746 (1994); F. Mertens and R. Imbihl, Nature 
370, 124 (1994). 

[4] H. Meinhardt, Models of Biological Pattern Formation 
(Academic Press, London, 1982); J. D. Murray, Mathe- 
matical Biology (Springer, Berlin, 1990); A. Koch and H. 
Meinhardt, Rev. Mod. Phys. 66, 1481 (1994); H. Mein- 
hardt, Nature 376, 722 (1995). 

[5] V. Castets, E. Dulos, J. Boissonade and P. De Kepper, 
Phys. Rev. Lett. 64, 2953 (1990); Q. Ouyang and H. 
L. Swinney, Nature 352, 610 (1991); K. J. Lee, W. D. 
McCormick, H. L. Swinney and J. E. Pearson, Nature 
369, 215 (1994); L R. Epstein and V. K. Vanag, CHAOS 
15, 047510 (2005). 

[6] K. Kaneko, Prog. Theo. Phys. 72, 480 (1984); J. P. 
Crutchfield, Physica DIO, 229 (1984); T. Yamada and H. 
Fujisaka, Prog. Theor. Phys. 72, 885 (1984); K. Kaneko, 
CHAOS 2, 279 (1992). 

[7] J. D. Keeler and J. D. Farmer, Physica 23D, 413 (1986). 

[8] L Waller and R. Kapral, Phys. Rev A 30, 2047 (1984). 

[9] V. A. A. Jansen and A. L. Lloyd, J. Math. Biol. 41, 232 
(2000). 

[10] E. A. Jackson, Perspectives in Nonlinear Dynamics, Vol. 

2 (Cambridge University Press, 1990). 
[11] P. Rohani, R. M. May and M. P . Hassell, J. Theor. Biol. 

181, 97 (1996). 
[12] R. Kapral, Phys. Rev. A 31, 3868 (1985). 
[13] T. Yamada and H. Fujisaka, Prog. Theor. Phys. 70, 1240 

(1983). 

[14] H. Chate and P. Manneville, Physica D 32, 409 (1988). 
[15] H. Chate and P. Manneville, Europhys. Lett. 17, 291 
(1992). 



[16] A. Politi, R. Livi, G. L. Oppo and R. Kapral, Europhys. 

Lett. 22, 571 (1993). 
[17] M. Yokozawa, Y. Kubota and T. Hara, Ecol. Model. 106, 

1 (1998); M. Yokozawa and T. Hara, Ecol. Model. 118, 

61 (1999). 

[18] M. Doebeh and T. Killingback, Theor. Popul. Biol. 64, 
397 (2003). 

[19] H. Meinhardt and A. Gierer, J. Ceh. Sci. 15, 321 (1974); 

J. Collier, N. Monk, P. Maini and J. Lewis, J. Theor. 

Biol. 183, 429 (1996). 
[20] J. R. Sanchez, J. Gonzalez-Estenez, R. Lopez-Ruiz and 

M. G. Cosenza, Euro.Phys.J. ST 143, 241 (2007). 
[21] J. Gonzalez-Estenez, M. G. Cosenza, R. Lopez-Ruiz and 

J. R. Sanchez, Physica A 387, 4637 (2008). 
[22] J. Gonzalez-Estenez, M. G. Cosenza, O. Alvarez- Llamoza 

and R. Lopez- Ruiz, Physica A 388, 3521 (2009). 
[23] J. Hofbauer, V. Hutson and W. Jansen, J. Math Biol. 

25, 553 (1987). 

[24] J. M. Gushing, S. Levarge, N. Chitnis and S. M. Henson, 

J. Diff. Eqns. AppL, 10, 1139 (2004). 
[25] R. M. May, Nature 261, 459 (1976); M. J. Feigenbaum, 

J. Stat. Phys. 19, 25 (1978); M. J. Feigenbaum, J. Stat. 

Phys. 21, 669 (1979). 
[26] P. J. Davis, Circulant Matrices (Chelsea Publishing, New 

York, 1994). 

[27] J.-F. Sadoc and R, Mosseri, Geometrical Frustration 
(Cambridge University Press, Cambridge, 1999); D. Nel- 
son, Defects and Geometry in Condensed Matter Physics 
(Cambridge University Press, Cambridge, 2002). 

[28] D. Chowdhury, Spin Glasses and Other Frustrated Sys- 
tems (Princeton University Press, 1987); J. Mydosh, Spin 
Glasses (Taylor & Francis, 1995). 



